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Abstract 

We consider an air bubble in water under conditions of single bubble sonoluminescence (SBSL) 
and evaluate the emitted sound field nonperturbatively for subsonic gas-liquid interface motion. 
Sound emission being the dominant damping mechanism, we also implement the nonperturbative 
sound damping in the Rayleigh-Plesset equation for the interface motion. We evaluate numerically 
the sound pulse emitted during bubble collapse and compare the nonperturbative and perturbative 
results, showing that the usual perturbative description leads to an overestimate of the maximal 
surface velocity and maximal sound pressure. The radius vs. time relation for a full SBSL cycle 
remains deceptively unaffected. 
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1 Introduction 

This extensive study of the sound produced in single bubble cavitation is prompted by the recent 
interest in sonoluminescence (SL), i.e. the conversion of externally applied sound in a liquid into 
light [^. It has been generally accepted that the SL light emission is intimately connected to the 
dynamics of the gas-liquid interface, though the specific mechanism producing the light flashes has 
not been uniquely identified. We present here a nonperturbative study of the sound radiated from the 
subsonically moving phase boundary between the gas bubble and the liquid in which it is immersed. 
In our approach we rely on the standard treatise of the subject [^. Apart from the importance of 
sound emission as the primary damping mechanism under the extreme conditions of sonoluminescence, 
the sound signal can be used as another diagnostic tool to study the bubble dynamics. For detailed 
accounts of related previous theoretical studies of the cavitation sound see Refs. i, @, including 
references to earlier work. We proceed here in two distinct steps: 

• To begin with, the motion of the gas-liquid phase boundary is assumed to be given and we 
determine the resulting sound field. 



Secondly, we consider the feedback effect of the radiated sound field onto the phase boundary 
dynamics, i.e. we describe the motion of the bubble selfconsistently. 



The (multi-)bubble cavitation has been studied intensely for some time, because of important 
applications and interesting underlying fluid dynamics 0. The discovery of well controlled 
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experimental conditions by Gaitan and Crum permitting the study of a single gas bubble over long 
periods of time, together with the ability to drive it externally by acoustic waves, has focused the 
experimental interest on the single bubble sonoluminescence (SBSL) phenomenon. The remarkable 
finding here is that light is produced in very short pulses Q: a million photons of several eV energy are 
emitted within 50 ps by apparently about 10^'^ atoms/molecules from a gas bubble of sub-/xm radius. 

The standard theoretical tool in the study of the bubble surface dynamics is the Rayleigh-Plesset 
(RP) equation, representing the motion of the gas-liquid interface by the time dependent (spherical) 
bubble radius R{t) . Several mechanisms which damp the collective motion of the bubble-liquid system 
have been incorporated in modern approaches: viscosity of the liquid and sound damping Effects 
of the 'rectified' mass diffusion |||, ^ and heat conductivity, see e.g. Ref. |1C], play also an important 
role in establishing the parameter regime for stationary SBSL, but will not be further discussed here. 
In the SBSL cycle the sound emission produced during the violent bubble collapse turns out to be the 
most important energy loss mechanism [Q, Less than half of the energy in the bubble is dissipated 
by the viscous friction. Only a comparatively tiny fraction of the energy delivered to the bubble by 
the external sound wave is radiated in form of the visible light H, |l^ . 

We will use the linear acoustic approximation (see Section 2.1), which describes the subsonic propa- 
gation of density perturbations in a fluid neglecting the effects of sound wave dispersion and absorption 
Q. This approximation does not allow us to study situations involving shock wave formation in the 
liquid. In Section 2.2 we solve the sound wave equation for the velocity potential exactly for any given 
R{t) . Our derivation is formally exact provided the assumptions of non-dissipative fluid dynamics 
are valid. Our expressions account in full for retardation and allow an arbitrary incoming sound field. 
The key result is the exact form of the pressure amplitude radiated from the moving bubble surface 
into the surrounding liquid, see Eq. (p8|). In order to compare with earlier analytical and numerical 
work we also derive in Section 2.3 the perturbative expansion of our general results in powers of R/c, 
i.e. bubble wall velocity divided by liquid sound velocity. 

In the subsequent Section 3 we describe a model of the bubble interior which is motivated by the 



approximately homologous dynamics observed in numerical hydrodynamic simulations |12]: the gas 
density distribution and velocity field have been found to show a simple scaling behavior, in terms 
of the scaled radius variable ^ = r/R{t) and with v{r,t) ~ ^R, respectively. These numerical results 
have recently motivated development of a semianalytic approach for the coupled bubble-liquid system, 
based on a Lagrangian variational principle ||l^, 0], of which a simpler variant will be used here. This 
allows to assess the order of magnitude of the feedback effects of the bubble interior dynamics. 

With this preparation we can proceed to the second step indicated above in Section 4, i.e. the 
derivation of a selfconsistent bubble equation of motion. There exists already a variety of equations 
generalizing the Rayleigh-Plesset equation, see e.g. [^, ^, |^, and more recent ones which consider 
the case of the externally driven damped nonlinear bubble oscillations in the SBSL parameter regime 
m, In our selfconsistent approach (Section 4.1) we make extensive use of the results obtained in 
Section 2 . We incorporate the external sound field in Section 4.2 according to standard experimental 
conditions realized to date. In Section 4.3 we derive a nonperturbative equation, Eq. (|89|), which 
is formally valid to all orders in R/c (< 1), however, subject to the linear acoustic approximation 
scheme. We incorporate the bubble interior using the vanderWaals equation of state and allow also 
for a realistic equation of state (EoS) of the liquid. Our derivations are based on the Navier-Stokes 
equation and, thus, take the effects of the liquid and gas viscosities consistently into account. 

We illustrate our formal results in Section 5 by studying numerically the properties of an air 
bubble in water externally driven by an ultrasound field with system parameters in the typical SBSL 
regime ||l], |ll| . We show in detail there that the various (nonperturbative) corrections incorporated in 
the bubble equation of motion generally result in sizeable corrections of the maximal bubble surface 
velocity during the collapse. They tend to make the bubble collapse less violent than seen in previous 
perturbative descriptions. In Section 5.3 we compute the outgoing compression wave emitted from 



The Sound of Sonoluminescence 



3 



the bubble surface into the surrounding hquid and consider the vahdity of the acoustic approximation. 
We present final remarks and conclusions accompanied by a brief summary of our work in Section 6 . 

2 Sound Emission from Phase Boundaries 

We reexamine here the sound radiation originating from a moving phase boundary between two non- 
viscous fluids, separated by an idealized wall of vanishing thickness. References to earlier work on this 
subject are surveyed in Refs. [§, |6|, 0]. Related issues were studied recently in Refs. H], §]. In distinction 
to these approaches, we derive exact expressions for the sound field emitted from a spherical boundary 
undergoing arbitrary motion. We begin in the next subsection by recapitulating some aspects of 
nonviscous hydrodynamics which lead to the (acoustic) approximation scheme used here in describing 
the sound generation and propagation. 

2.1 Acoustic Approximation 

The linear acoustic approximation is based on the assumption that the sound field causes only small 
perturbations of the ambient state of the fiuid characterized by a vanishing velocity field and constant 
density and pressure Ref. ||2|. Thus we shall not consider here nonlinear phenomena, such as the 
dispersion and absorption of sound waves in the medium and, in particular, the formation of supersonic 
shock waves. We shall see quantitatively in Section 5 that this scheme is fairly well satisfied, though 
the amplitude of the generated sound field typically reaches up to 10^ atm near to a SBSL bubble. 
The basic equations of nonviscous hydrodynamics are: 

1. The continuity equation for the fluid density p , which assures the conservation of mass: 

dtp + v- (pv) = , (1) 

where v denotes the local fluid velocity. 

2. The Euler equation, arising from Newton's equation of motion for a small fluid cell: 

pDtv = - VP , (2) 

where P denotes the local pressure in the nonviscous fluid and the comoving derivative is deflned 
hy Dt = dt + V ■ V . Later we will consider the Navier-Stokes generalization for viscous flow, see 
Eq.(|D. 

3. In the absence of entropy production the set of equations is closed by the EoS relating the 
pressure to the density. For example, the ideal gas EoS is: 

P[p] =Kp^ , 7 ^ Cp/Cv , (3) 

where the adiabatic index 7 denotes the ratio of the specific heat at constant pressure and 
volume, respectively. 

Had we allowed that the sound waves disturb the medium and generate entropy, we would require an 
EoS with two variables, e.g relating the pressure to the mass and energy densities. In this case we 
would have to introduce additionally the energy conservation in differential form, in order to close the 
set of equations. 

In the following we study small-amplitude, acoustic perturbations of an ambient state of the fluid 
(liquid, gas) characterized by the ambient solutions to Eqs. (||)-(|3|), Pq; vq, pq. In order to obtain 
equations governing the perturbations P* , v*, p* , we set: 



P = P0+P* , v = Vo + v* , p = pQ+p* 



(4) 
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Inserting these expressions into Eqs. (|Tl)-@, a coupled set of equations results which relates various 
powers of the perturbations. The description of perturbations can be further simplified by assuming 
a particularly simple ambient state: 

Po{x,t) = Const. , vo{x,t) = , po{x,t) = Const. , (5) 

i.e. a homogeneous state independent of time (quiescent state). Then, the equations linearized in the 
perturbations become: 

dtp* + poV-v* = , (6) 

podtif* = -VP* , (7) 

P* = jKpl'p* = c'p* , c'^^ , (8) 

where we made use of a Taylor expansion in order to obtain Eq. and introduced the (constant) 
sound velocity c for the ambient fluid state. Combining Eqs. (|6|)-(|8), we obtain the well known wave 
equations: 

d'^P* - c^AP* = , d'^p* - c^Ap* = . (9) 
Next, we introduce the velocity potential cp . Consider the curl of Eq. (^): 

podtV XV* = -V X VP* = , (10) 

i.e. V X w *(x, t) = C{x) . Thus, if the vorticity of the velocity field perturbation v* vanishes initially, 
then it remains zero always. In this case we can write: 

^* = V0 , p* = -podtci) , (11) 

where the relation between the velocity potential and the pressure perturbation is chosen such that 
the linearized Euler equation, Eq. (0), is automatically satisfied. Inserting the EoS, Eq. (^), in the 
linearized continuity equation, Eq. (P), and using Eq. (11), one obtains the wave equation for the 
velocity potential: 

d!^-c^A^ = . (12) 
For simplicity, we suppress from now on the index * introduced in Eq. (^. 

2.2 Outgoing Spherical Wave Dynamics 

Assuming spherical symmetry and using spherical coordinates, the velocity potential is seen to satisfy 
the usual scalar spherical wave equation: 



r2 

= r{dt(l) 52(r0)} = d?(r(l)) - c^d^lrcj)) . (13) 

r 

The generic solutions (regular at infinity) are out- and ingoing spherical waves: 

(l)out=r-^f{t-r/c) , ^in = r-^git + r/c) . (14) 

For the spherically symmetric bubble wall located at the position r = R{t) the boundary condition is: 

v{R, t) = dr<P{r, t)\,=R^t) = R{t) , (15) 

where we introduced the notation X = dX/dt . In general, the solution of the wave equation, Eq. (14), 
will be composed of a linear superposition of in- and outgoing spherical waves. We shall allow for the 
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possibility of an incoming acoustical wave, <pin: as part of the external driving sound field, and 
rewrite Eq. (^) in more detail: 

11- 

dr4'out\r=R = -^fi^ " ^A) " ^/'(^ ~ ^/^"l = R- dr(pin\r=R = Ti , (16) 



where (l)out has been substituted by an outgoing spherical wave, cf. Eqs. (|14D ; henceforth we use the 
abbreviation f'{x) = df/dx. The external field (pin is kept arbitrary in the derivations to follow, but 
will be incorporated according to a physically relevant example in Section 4.2 . 

In principle, there will be also a sound field generated by the moving phase boundary which will 
travel into the interior, to be reflected at the center of the bubble (r = 0), and then return. Its 
description, i.e. pertaining to the bubble interior, requires a suitable modification of the present 
approach, see Section 3 . 

With TZ{t) given by other dynamics, our objective now is to determine implicitly the outgoing 
velocity potential (pout ■ For this purpose we recast Eq. ( p!6[) into the form of an ordinary linear first- 
order differential equation with known time dependent coefficients: 

1 1 1 ■ ■ 

-^F ^F = n , (17) 

cRl-R/c ^ ^ 

where we used of the following substitutions and simple relations: 

t-{t) = t- R{t)/c =^ i- = l-R{t)/c , (18) 

df dt^l f{t^) 

f(t-Rc) = — — —^ = - -r— , 19 

dt- dt t- 1 - R/c 

F{t)^fiUt)) =^ F{t) = f{Ut)) . (20) 



In Eq. (H) we note the retardation factor originating in Eq. (p!9[). The general solution of Eq. (|l7|) is 
elementary: 

Fit) = Fo^^e-^^o^'^ 

-R{t) f' dt"[c- R{t")]TZ{t")e'-^^"'^''T^ , (21) 
Jto 

where Fq = F{tQ) denotes the integration constant. It is fixed by the requirement that no outgoing 
wave should be excited before the bubble wall is set into motion at Iq . Therefore, with TZ(t < to) = 0, 
we find Fq = . The resulting constraint t > to is most conveniently implemented by a step function 
factor 6{t — to) . 

Using Eqs. (p!4[), (20), we can recover from the above solution the outgoing velocity potential. To 



do this we will need the inverse function t(t-) = t of t-{t), Eq. (18). Note that t is a single- valued 
function, if and only if t-{t) = t — R{t)/c is either strictly increasing or strictly decreasing with time. 
For subsonic bubble wall motion t_(t) is strictly increasing, since 1 — R/c > 0, and t is strictly 
increasing in this case. This allows us to perform variable substitution at will and we obtain: 

cPoutir,t) = -F{i{t-r/c)) = - to)^ T (it"[c - i?(t")MO e" , (22) 

r r Jto 

s = i{t-r/c)=t-^^^^^ . (23) 



6 



H.-Th. Elze, T. Kodama and J. Rafelski 



The retardation effect seen here can be evaluated more exphcitly. For example: 

iit-r/c) « ^_ r-R{t-r/c) ^ ^ ^ j^^^^ ^ ^24) 

i{t-r/c) ^ , r^R{t) . (25) 

c 

In order to calculate the velocity and pressure fields, we require the derivative of t: 

t{t - R{t)/c) ^ t =^ i'{t^) = ■ 1 . (26) 

1 — R{t{t-))/c 



Then, differentiating Eq. (p^), we obtain (s > to): 

v{r, t) = dr[(pout{r, t) + 4)in{r, t)] 



+ +dAin{r,t) , (27) 



r R r r 



i{t~r/c) 



P{r, t) = -pLdt[(l)out{r, t) + 4>in{r, t)] 

= ^\Rn+^} -pMin{r,t) , (28) 

where all functions within the braces are to be evaluated at the indicated time-argument t{t — r/c) = s; 
henceforth pi denotes the ambient density of the liquid. We recall that P is the pressure perturbation 
due to the generated sound field, to be added to the ambient pressure. 

2.3 Series Expansion for the Sound Field 

The results obtained for the velocity potential and the velocity and pressure fields, Eqs. (21)-(p8|) in 
Section 2.2 , involve an integration over the history of the bubble wall motion. Thus the corresponding 
nonlinear differential equation describing it selfconsistently (Section 4) could be expected to become 
nonlocal in time (see, however, Section 4.3). In any case, for sufficiently slow motion a local approx- 
imation can be justified, see Refs. ^, and earlier references therein; this approach leads to a 
popular generalization of the RP equation. Therefore, our next objective is a systematic expansion in 
powers of R/c of our general expressions, in particular of the velocity potential, Eq. (|2^). 
Let us consider an integral of the form (in our case j = c/ R , / , h oc {1 — R/c)TZ): 

I{s)= r dt"h{t")exp{- r dt'j{t')} , (29) 
Jto Jt" 

with j > and s >tQ , cf. Eq. (21). For a sufficiently large integrand in the exponential the integration 
over t" will be effectively limited to a small range next to the upper limit s . This suggests the following 
expansion: 

lis) = rdt"hif)exp{- rdt'{j{t')-j{smexp{jis){t"-s)} 
Jto Jt" 



n=0 



oo 



1 7>-{n)/ ^ d"" 1 - exp{i(s)(to - s)} 



^^nl ^ 'dj^ j{s) ^ ' 



The Sound of Sonoluminescence 



7 



where: 



hii)eM- I dt'ij{t')-j{s))} 



(31) 



t=s 



We remark that the exponentially small correction term (for s ^ Iq) on the r.h.s. of Eq. ( pO| ) can be 
attributed to transient contributions to the integral or, rather, to the function F, Eq. (^), which are 
due to the switching-on of the bubble wall motion considered in Section 2.2 . Since we will be mainly 
interested in the nonlinear oscillatory motion of the bubble wall when all transients have died out, we 
can neglect this term for sufficiently late times. Thus, we obtain: 



(32) 



n=0 



which presents the starting point of our approximate evaluation of the velocity potential. 

Employing Eqs. (|2l|), (|2^) from the previous section together with Eqs. (|29|), (|3^) above, we obtain: 



f>out{r,t) 



(33) 



E( T^) — (l-^W)exp{(5-t)ij77T+ / dt^} 



n=0 



dV' 



R{s) 



t=s 



7^ <^ 1 - 2-1 + + 2-^1 + 2^ + + n-rr^ + O ({R cf 

r \ 2 RTZ ^c^^ RTZ R^ 2 R^n 



where as before s = t{t — r/c) see Eq. (ps]). The expansion has to be carried out up to order n = 4 
in order to acquire in the R/c expansion all second order terms. Our result reduces to Eq. (18) of 
Ref. Q at the 0{R/c) obtained there, provided we replace i) TZ{t) — > R{t), i.e. we ignore effects of 
an incoming sound field, cf. Eq. ([l^), and ii) i{t — r/c) — > t — [r — R{t)]/c . The latter approximation 
amounts to neglecting already some 0{R/c) corrections in Eq. ( p^ ) due to retardation effects. 

To conclude this subsection, we calculate the velocity and pressure fields generated by a driven 
oscillating bubble wall, i.e. the sound emitted, by combining Eqs. (^), (^) and ([2^), respectively, 
and Eq. (H): 



v{r,t) = dr(pin{r,t) + 



Rl'n 



(34) 



+ 2- 



RRK 




1 RTZ 
2RR 



c RTZ R^ 2 R^TZ ^ ' ' 



which illustrates nicely that our expansion reproduces the incompressible fluid (c oo) limit; fur- 
thermore: 



P(r,i) 



+ ^RRTZn + l^ 
r [ 2R'JZ 



R,^ RIZ RR IR^TZ, 
- 1 + 2^ + ^ + -^-^ + 
c ' RTZ R^ 2 R^TZ 




(35) 



where the inverse of the retarded time function has to be inserted everywhere as indicated before. 
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3 Homologous Bubble Interior Dynamics 

In section 4 we want to derive the equation of motion of the phase boundary (bubble surface) self- 
consistently. For this purpose we need to understand the dynamics of the interior of the bubble. 
In distinction to Section 2 , presently we have to take into account that the density, pressure, and 
velocity fields inside the bubble may change by several orders of magnitude during different phases 
of the (periodic) bubble motion. Therefore, a linearization of the hydrodynamic equations of motion 
around a homogeneous and time independent ambient state, as performed in Section 2.1 , is not 
applicable here. 

In the hydrodynamic studies of the interior motion the emphasis has been to understand the 
development of the extreme conditions inside the bubble ||l5| , |l^. Chu noted that for an important 
part of the SBSL cycle the bubble motion is (nearly) homologous, i.e. the shape of the density 
distribution and the velocity field of the gas scale in an appropriate way with the radius R and surface 
velocity R ||l^. The numerical simulations in the SBSL parameter regime indicate a homologous 
contraction (and its stability) until the onset of shock wave formation and independently of the details 
of the EoS used. This stability of the homologous motion may be related to the fact that only for the 
homologous motion there is no energy loss due to viscosity inside the bubble (see Eq. ( |60| ) and below). 

We solve the continuity and Euler equations, Eqs. (|l|) and (§) of Section 2.1 respectively, by ap- 
proximating the density and velocity fields with suitable scaling functions. We neglect dissipative 
effects in the bubble, such as viscosity, heat conduction, and radiation transfer especially. Presently, 
we are interested in a semiquantitative analysis of the bubble interior. This seems sufficient for our 
later derivation of the overall bubble dynamics: the motion is largely determined by the dynamics of 
the fluid, since only a minor fraction of the relevant total energy resides in the interior at any time. For 
a detailed understanding of the microscopic processes leading to sonoluminescence inside the bubble, 
however, the dissipative effects will be essential ingredients. 



It turns out, see Refs. |13, 14 1, that for a time independent shape of a properly scaling density 



distribution the velocity potential is: 



^ir,t) = ^|r2 , r<R , (36) 
in terms of the bubble radius R{t), which yields: 

v = dr^=^r , M = \(^^-{^fy' , (37) 

i.e. a linear radial velocity profile. The continuity equation, Eq. (|l]), takes the form: 

= dtp+ ^drir'^pv) =(^Q + 3^0^ d + Qdtd\i: , (38) 

where we introduced the homologous ansatz: 

p{r,t) = Q{t)d{i,t) , (39) 



with ^ = r/R{t) denoting the scaling variable. Here we find indeed that the velocity potential, Eq. (|36D, 
is consistent with a time independent density profile function, i.e. 

d{i:t)=d{i) , (40) 



and the overall density factor of Eq. (39), normalized to N particles inside the bubble: 



N 

Q{t) = . . (41) 

^^R^{t)!^died{i) 
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We choose d{l) = 1 in the foUowing, i.e. p{R,t) = g{t) . We refer the reader to Refs. ||T^, 14] for the 
variational improvement of the strictly homologous dynamics (see and references therein) which 
we will pursue here, with a time independent yet scaling density profile function. 

In order to determine the scaling function d(^), we consider the integral of the Euler equation, cf. 
Eq. (0), which for our radially symmetric case yields: 



r I 

dr -drP . 
R P 



In the absence of entropy S production the integral on the r.h.s. is the enthalpy: 

1 rR dP 

P P 



rR I i-R dP P 

/ dr-drP= / \s = h(R) - h(r) , h = £+- 

Jr P Jr P P 



(42) 



(43) 



where £ = E/N is the specific energy per particle. 

During the highest compression particle densities reached in the bubble interior are similar to those 
in the liquid outside. We thus use here for the gas an adiabatic EOS, cf. Eq. (^, with vanderWaals 
hard core corrections: 



Pip) = Po 



Pa - PO 
Pa- P 



Pn 



p/po 



'^-P/Pa 



(44) 



Here po (^ pa) is the density at which the bubble interior is at the ambient pressure Pq and p^^ = 
^o? /N in terms of the vanderWaals excluded volume. The adiabatic index 7 for ideal monatomic 
(diatomic) gases is 5/3 (7/5) . Related EoS have been applied to calculate the gas pressure before 



Performing the integral in Eq. (|4^ ) using Eq. (44) and employing the scaling form for the velocity 
potential, Eq. (^), we obtain = r / R{t)): 



1) 



Pn 



p/po 

po \1- p/pa 



Po \Po 
Pq ( Pa 



7-1 



7-1 
7 

7-1 
1 



1 



1 - p/pa 

+ o(p/a 



+ 



1 



7-1 



R 



R 



Pa \PoJ (1 - p/paV 



[l + 0{l-p/pa)) 



(45) 
(46) 

(47) 



R 



where Eq. ( ^6| ) and Eq. (|47|), respectively, represent the low- and high-density limits of the r.h.s. of 
Eq. (H). 

Obviously, the vanderWaals correction introduces an additional scale which makes it difficult to 
find a universal density shape function. Therefore, we consider the low- and high-density cases in 
turn: 



Low density. Inserting Eqs. (|39|)"(p|) into Eq. (|4^), we obtain: 



ipoRR ( PO 

Q{t)) 



27 P( 



(48) 



Since the l.h.s. here only depends on t and the r.h.s. only on ^, we have introduced a constant 
parameter q\ characterizing the time independent density profile. Solving the r.h.s., we find: 



di{0 = (l + qx{l-e) 



1 

7-1 



(49) 
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We note that the l.h.s. of Eq. (^) determines the value of the parameter qi, and thus we can 
verify if indeed it is time independent. Prom this we conclude that this low-density (p <C pa) 
functional form of the density profile function is approximately valid for inertial motion with R ~ 
or, more generally, when the bubble wall accelerates (decelerates) proportional to the internal 
pressure acting on it, P{R,t) « Po{Q{i) / PoV ■, that is whenever the ratio {gRR}{t)/ P{R,t) is 
approximately constant. 

• High density. In this limit, inserting Eqs. (|3^)-(^l|) into Eq. (|47|), we obtain: 

1 fpoY PaRR il-Qit)/par+^ _ _m-l 

where again the l.h.s. is only a function of time, while the r.h.s. is only a function of ^, and thus 
we can introduce the constant g^. In deriving Eq. (|50D, we employed an additional expansion for 
d(^) fa 1 on the r.h.s. which is consistent with the high-density case (p ~ pa) under consideration. 
Solving the r.h.s. for the density profile function, we obtain: 

dhiO = ^ + Qh{l - f) , (51) 

while the l.h.s. determines the value of the parameter ^-s before. The limit of validity of this 
result is controlled by the requirement that the expression on the l.h.s. remains time independent. 
Note that the high- and low-density profiles agree with each other for 7 = 2; however, 7 = 3/2 
and 5/2 for mono- and diatomic gases, respectively. 

We remark that the numerical values of qy and q^, in general, vary depending on the dynamical regime, 
and always ^ > — 1. It is important to realize that Eqs. (^8|) and ( ^0|) present implicit equations for q\ 
and q\i, respectively. The reason for this is that g{t) because of the normalization in Eq. ( p| ) depends 
on the integrated density profile and thus on the respective q by Eqs. (^9D, (|5l|). In practice, these 
'constants' will be treated as adiabatically changing parameters to be computed selfconsistently from 
the equations derived here. 

We observe that as the bubble goes through the cycle determined by the driving pressure, the 
shape of the matter distribution inside the bubble changes, controlled by q, which remains nearly 
constant during much of the cycle and changes rapidly in sign near to the bubble collapse/bounce, 
whenever R = 0, where the l.h.s. of Eqs. (|4^), (^8|), or (^) vanishes, and q = 0, assuring in this 
transition instant a (nearly) homogeneous density distribution. For negative values of q the highest 
density is found at the surface, corresponding to a collapsing phase; similarly, positive values of q 
with the highest density at the center correspond to an expanding phase of the bubble motion. We 
found that there is an approximate fixed point r/R = in the motion, where the average density 
corresponding to a homogeneous bubble is maintained during the entire cycle: 

p^ = p{C^R,t)c^N/{4TT/3R{tf , C*« 0.78. (52) 

This will turn out to be a very useful property which we shall exploit deriving specific numerical 
results in Section 5.2 , where we compare the bubble dynamics for a homogeneous and a homologous 
interior, respectively. 

4 Derivation of the Bubble Equation of Motion 

Having studied the effects of a given bubble wall motion on the exterior liquid and the bubble interior 
in Sections 2 and 3 , respectively, we now turn to the question how this motion can be determined 
selfconsistently. That is, assuming an external driving sound field (pin, we want to derive an equation of 



The Sound of Sonoluminescence 



11 



motion for the bubble radius i? as a function of time. We also take the bubble interior into account. Its 
dynamics may have little influence on most of the cycle of stable bubble oscillations. However, in order 
to improve the understanding of the final stage of the violent collapse [13, 14], it has to be considered. 



Besides the sound emission from the moving phase boundary, which constitutes a major energy loss 
mechanism, the effect of viscous damping will also be incorporated here, which is of a similar order 
of magnitude as acoustic damping. Despite the fact that our considerations so far were based on the 
assumption of nonviscous fluid dynamics, i.e. approximately free wave propagation in particular, we 
can include dissipative and driving forces into the equation of motion for R{t) by deriving it from 
the Navier-Stokes (momentum balance) equation. Thus, only the backreaction of viscosity (and other 
dissipative transport effects) on the wave propagation and velocity potential cj) is neglected here. 

4.1 The (Generalized) Rayleigh-Plesset Equation 

We begin with the Navier-Stokes equation generalizing the Euler equation, see Eq. (§), to the case of 
viscous fluids For our spherically symmetric situation this equation can be written in the form: 

dr [dtcP + ^{drcPf^ = -UrP + ( ^^^'(r^)) , (53) 

where (/>(r, t) denotes the appropriate velocity potential, i.e. for the interior or exterior region, and 
r/ = |r7<j + r/fe is the relevant combination of shear and 5ulk viscosities of the fluid entering here. 
Similarly, the pressure P has to be specified differently according to whether the interior or exterior 
of the bubble is considered. 

Due to the phase change at the bubble surface, there arises a nontrivial boundary condition relating 
the normal components of the stress tensor S inside (Gas) and outside (Liquid): 

n-T.G-n = n-T.L-n-'^ , (54) 

where the term oc a denotes the 'pressure' contribution due to the (liquid) surface tension for the 
liquid/gas interface and the viscous stress tensor has to be evaluated for the respective fluid under 
consideration 0: 

n ■ ■ n = \ —P + rjdrV — —v\ , (55) 
where f) = ^r]s — ^rjb, rj was defined after Eq. (|5^), and v{r,t) = dr4>{r,t), as usual. Then, employing 



the boundary condition (15) and the velocity potential ( pq ) for the gas phase, we obtain instead of 
Eq. (^) more explicitly: 

- Pg{R) + iVG - %)| = -PUR) + VLdrV\r=R " ^^f " ^ • (56) 



Using Eq. (|26|) and the velocity field (|34D , we calculate drv\ii for the outer liquid: 



drV{r,t)\r=R = -2^-R-U[-2^ + -dt + {-fdt]dr^,n{r,t) 

H \ H c c 



r=R 



- 2^i^f f 1 + 2^ + a + O UR/cf) , (57) 

where the first term on the RHS corresponds to the incompressible fluid limit (c oo) in the absence 
of an external driving sound field. 
Several remarks are in order here: 
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The bubble gas pressure at the surface, Pq {R) = Pg {pg{Rt t)), is determined by the van der Waals 
EoS Eq. (|^) together with pg{R, t) = Qcit), according to Eqs. (|39|), ( |4l| ) for the homologous den- 
sity function. The liquid vapor pressure inside the bubble is not considered at present. 

We observe that rjc — fjc = 3r]i,G, in terms of the bulk viscosity of the gas, which is generally 
very small compared to liquid viscosities. The shear viscosity does not contribute here at all, 
cf. also the remark after Eq. (|6^) below. For completeness, we allow for viscous damping in the 
gas inside the bubble, even though it might at best become important during the rapid collapse 
and high compression phase. 



We recall that by definition of the inverse of the retarded time function, Eq. (pq), we have 
i{t — R{t) /c) = t ; therefore, the r.h.s. of Eq. (|57| ) has to be evaluated at the time t . 



In the following the boundary condition ( pq ) together with Eq. ( |57| ) will be employed to eliminate 
Pl{R) in terms of the other quantities which by now are explicitly calculated functions of R, R, R 
etc. In order to proceed, we convert the Navier-Stokes equation ( [5^ ) together with the boundary 
condition (]5^ ) into an ordinary differential equation for R{t). This can be achieved by integrating 
over the radial coordinate, i.e. the approach advocated earlier in Ref. for example. Symbolically, 
in obvious correspondence with the terms on the left- and right-hand sides of Eq. (|53D, we obtain: 



rR roo rR roo rR rOD 

{ + / }drV={ + / }drV + { + / }drV , 
Jo JR Jo Jr Jo Jr 



(58) 



where we split the integrations at r = R because of the different phases inside and outside. Making 
use of the results of the previous Sections 2 and 3, we proceed to evaluate each term of the generalized 
Bernoulli equation (|58| ) in turn. 

Using the velocity potential (p6|) for the bubble interior, we obtain: 



^ drV = drdr \dt(t> + \{dr(ff^ = \rR , (59) 

f^drV = r dr'^dr(-d'^(r(j))] =0 . (60) 
Jo Jo P \r ) 



We observe that the homologous scaling solution for the bubble interior, </> oc r^, does not give a 
contribution to the integrated viscous force term here. This remarkable result suggests that viscous 
forces drive the bubble motion to the homologous limit which we discussed in Section 3 and which is 
in this respect quite unique |Q, 14 1. 

Furthermore, employing the van der Waals EoS (^), we calculated the enthalpy integral for the 
gas inside the bubble before, cf. Eqs. (43), (|45|): 



/ drV = -f dr -drPG = - Pg{pg) , 
Jo Jo p (7 - 1)/0G 



where we rewrote our previous result of Eq. (|45|) in terms of the van der Waals expression for the 
pressure and denote by he the gas enthalpy. 

Next, we turn to the evaluation of the corresponding integrals for the exterior region. The first 
integral follows immediately from our derivations in Section 2 : 

foo poo r 1 1 1 1 

/ drV= drdr\ dtcl) + -{dr(t>f = —P(.R, t) - -R^ , (62) 
Jr Jr V 2. } pi 2 



where P(r, t) denotes the pressure field calculated in Eq. (p^) and evaluated approximately in Eq. (|35[); 
henceforth we denote by pL the ambient density of the liquid. The second term on the r.h.s. of 



The Sound of Sonoluminescence 



13 



Eq. (^) simply follows from the boundary condition (|l5| ) and the assumption i;^ ^ (sufficiently fast) 
for r — > oo . 

In order to calculate the enthalpy integral for the liquid, we employ an EoS which gives a realistic 
description for many liquids (5|: 

PLip) = (Po + Pi) [yj ^ ^ = ' (63) 

where Pq denotes the ambient pressure and n, Pi are parameters depending on the liquid (e.g. for 
water n = 7, Pi = 3kbar). Then, it is straightforward to obtain: 



oo 



a. v.- 1 *ia,P, = ^^^ (A^l m 



R Jb P n - 1 PL I V Po + Pi 

using P(r — > oo) = Pq , i.e. the ambient pressure, and where Hl denotes the liquid enthalpy. As 
mentioned before, the liquid pressure on the bubble surface, Pl{R), will be eliminated via the pressure 
boundary condition (p^). 

The last (viscous) term from Eq. (|5^ ) can be evaluated employing the approximation p ^ pL- 



r drV ^ r dr "^dr (-df.{rcP)) = -!II^-d^{rcl))\r=R = {2^ + a.t;|.=ij| , (65) 
JR JR p \r J pir PL [ P J 

i.e. neglecting a cross term between viscosity and compressibility of the liquid, which will be justified 
by the numerical results presented in Section 5.3 . 

It is worth while to recall that in Eqs. (^)-(p5|) all the functions R, R, R, etc., which appear on 
the r.h.s. , are to be evaluated simply at the time t . 

Collecting the essential results from Eqs. (|59[)~^ (|6^), we obtain the equation of motion for the 
bubble radius R(t) , i.e. /q°° drV = dr{V + V) , in the form: 

\RR-lR^ + ]_Pi^R^t) = hG[PG]+hL[PL]- — U^+drV\r=lA . (66) 
2 2 PL PL [ R J 

This equation is valid for any subsonic bubble wall velocity, it is nonperturbative in R/c. In fact, 
inserting here our general results of Section 2.2 for the sound field generated by the bubble, in particular 
v{r,t) and P(r, t) calculated in Eqs. (|27| ) and (|2^ ) respectively, we obtain an important generalization 
of the Rayleigh-Plesset equation. 

In order to illustrate the contents of Eq. (^), we proceed with a perturbative evaluation of terms 
involving the generated sound field making use of our results of Section 2.3 . Using in Eq. (|6^) the 
expansions in R/c, especially Eqs. (|35|), (57), we obtain the perturbative bubble equation of motion: 



-RR - + —P^ + 2P7t(l - -) + P:^(1 - 4-) - -RTZR - -R^ 7^ 
2 2 PL c c c c 

= haiPG] + hL[PL] - 2^ {[^ - ^^*]^«(^'*)}^^^ + O ((^/^)') • (67) 

We observe that for TZ = 0, i.e. when the bubble wall velocity is identical to the incoming acoustical 
velocity field at r = R{t), all terms involving 1/c-corrections vanish. However, on the r.h.s. the viscous 
damping term survives, which involves only the incoming sound field. This is to be expected, since no 
outgoing sound wave is generated, in agreement with the nonperturbative results for the velocity and 
pressure fields, v = Va = dr4>in and P = Pa = —pLdt'pin in this case. 
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We arrive at a more useful form of the bubble equation of motion by eliminating R from Eq. (| 
This can be achieved at the same order in R/c by using the equation at leading order (i.e. for c — > oo 
to calculate R and then R : 



RR=-\ ^[r.h.s.loo] - -RR 



PL 



-Pa + 2RVa + mi) a + RVa 



(68) 



where the term [r.h.s.|oo] denotes the r.h.s. of Eq. (|67|) in the limit c — > cxo and Pa = Pa{R,t 



Va = Va{R,t) , Pa = -^Pa{R,t) , Va = -^Va{R,t) , etc, from now on. Reinserting this expression into 
Eq. (^), we finally obtain the equation: 

■■ In/ '2va\ 3-n/ AR-Va\ 



-2RVa -RVail 



c / 3 c 



1 + (1 



+0 ((i?A 



1 



hG{PG\^hi\Pj^\-2 



T]L_ 
PL 



. 1 



2c 



-Pa 



--R PL 



(69) 



where we separated out the underlined terms, in order to facilitate the discussion of our result; the 
term oc d/dt\c^oo is to be read as: take this limit of the following expression before evaluating the 
derivative. Equation ( |69[ ) generalizes previously considered variants of the Rayleigh-Plesset equation. 
In comparison to earlier related work, particularly Refs. |^, ^, 13|, we draw attention to the following: 



The dynamics of the gas in the interior of the bubble as well as of the exterior liquid is fully 



incorporated in our derivation. It is based on the Navier-Stokes equation (53) supplemented 
by the boundary condition (56) relating the normal components of the stress tensor at the 
bubble surface. Neglecting the bubble interior amounts to setting the underlined numerical 
constants in Eq. (|6g|) to zero, equivalent to neglecting RR/2 in Eq. (p^), and setting hciPG] = in 
both equations. These apparently are sizeable corrections which were not considered previously. 
However, these terms cancel exactly, if the homologous dynamics of Section 3 is applicable for 
the bubble interior. Also, according to Eq. (|6l|), he = for a homogeneous bubble. 

The influence of the external acoustic driving field is consistently taken into account in our 
treatment of the sound field outside of the bubble. The terms oc Va, i.e. the incoming velocity 
field, or its derivatives, were neglected before; only terms involving the driving pressure Pa had 
been obtained. These additional terms are necessary to restore the boundary condition on the 
sound field at r = R{t), see the comment after Eq. 



has been employed here. For n ^ 1 the enthalpy 
- Po)/pl, which was used earlier instead. - The viscous 



The pressure boundary condition ( [Sq) contains terms which have been neglected before: 



A realistic EoS for the liquid, Eq 
in Eqs. dH), @) reduces to {Pl{R) 
term oc rji was not studied before. 



Pl{R) = PG(i?)-^-(3%G + 4r/.L)| 

+ dvsL + r]bL){^Va - -{dtVa{r,t)}B.) + O ({R/cf 
6 R c ^ 



(70) 
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where we used Eq. ( |57|) and inserted the shear and 6ulk viscosities for hquid and gas exphcitly. 
Except for a highly compressed or, rather, very hot gas Q, its bulk viscosity presents a negligible 
contribution in Eq. ([70|), being on the %-level as compared to liquid water at most. However, 
we also note the corrections due to the incoming sound field. 

If we discard altogether the modifications just discussed, we reproduce Eqs. (9) and (10) of Ref. 0]. In 
Section 5 we will explore numerically the effects of the corrections presently discussed. 



4.2 The Driving Sound Field 

Having further studies of single bubble sonoluminescence in mind, we determine here those quantities 
entering our results which depend on the external acoustic driving field. For typical experiments 
on single bubble sonoluminescence carried out to date the driving sound frequency is approximately 
1/ m 25 Khz, which translates into a wavelength of A ~ 5.9cm given the sound velocity in water of 
c = 1481 m/s. Relevant bubble radii are such that R/X ~ 10"'^ [Q, Therefore, we consider the 
external field in the long- wavelength limit. 

Following Ref. we assume a plane standing wave: (/> = A sm{k z) cos (ujt) , i.e. along the z-axis 
with Lo = 2ttv and k = 2tt / \ = uj /c . In the long- wavelength limit, with r ^ R <^ 1/k , the spherical 
component of the external field, determines the dominant driving force for spherically symmetric 
bubble oscillations. Expressing z = zq + rcosO, where zq denotes the presently irrelevant position of 
the center of the bubble, we can project out the spherically symmetric monopole component of the 
driving field to obtain: 

4'in{r,t) = — cos(u)t) / d(cos9) Po(cos9) sm(k[zQ + r cos9]) = Asm(kzo) cos{ujt) — ; , (71) 

2 7-1 kr 

with Pq = 1 denoting the appropriate Legendre polynomial. 

We recall that P{R, t) contains the contribution of the incoming sound field, i.e. the external 
"acoustic driving pressure" : 

Pa{r,t) = -pLdt(t)in{r,t) 

sin(A;r) n ^ ^\ /^-i ^ n \2 n,(n 



= PAsm{ut)^-^=PAs\Ti{u:t)[l--{krY + 0{{krf)'j , (72) 

cf. Eq. (p^). The amplitude of the acoustical pressure field is Pa = ujpLAsm.{kzQ) . Similarly, the 
velocity field v{r,t) has an external contribution: 

Vair,t) = dr(j)in{r,t) 

Pa , X f cos(kr) sm(kr)\ Pa . x /, . 



cf. Eq. (27), which also enters our results for the bubble equation of motion and pressure boundary 
condition. In general, both quantities. Pa and Va , are determined by the experimental set-up, in 
particular the external ultrasound sources which drive the bubble oscillations. 

Finally, we rewrite the relevant spherically symmetric part of the external velocity potential as a 
superposition of an outgoing and an ingoing spherical wave. We obtain: 

cPa 

0m(r, t) = ^-^{sin(a;[t + r/c]) - sin(a;[t - r/c])}/r , (74) 



which will be useful in the following section. 
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4.3 A Nonperturbative Bubble Equation of Motion 

We present in this section the bubble equation of motion resulting when we do not use the expansion 
in R/c which was employed in Section 4.1 , see Ref. for earlier related work. 

Since the velocity potential (f> obeys the spherical wave equation ([l^) for the presently considered 
case of a spherically symmetric gas bubble immersed in a liquid, we can immediately write down its 
most general form (see Section 2.2): 

(f)(r, t) = - if (tret) + g{tadv) " giUet)) , (75) 

r 

with arbitrary functions / and g and where: 

tret = t - r/c , tadv = t + r/c . (76) 

The terms involving g represent the external field, an example of which was considered in the previous 
section, see Eq. ([f^ in particular. Generally, the relative sign between the retarded and advanced 
contributions results from the requirement that the external field be regular at r = , i.e. the position 
of the center of the bubble. The term involving / presents the additional outgoing (sound) contribution 
to the velocity potential which is generated by the moving bubble wall. 

Then, our task is to determine (or eliminate) the functions / and g such that (p satisfies the velocity 
boundary condition (^): 

R = dr(P\R = -^{c<P + f'{tret)-g'{tadv)-g'{tret))R , (77) 

where f'{x) = df{x)/dx etc., and where Eqs. ([75|), ( [7^ ) were used to obtain its present form. Similarly, 
we rewrite the pressure boundary condition (|56|): 

Pl{R) - VL-d^r{r(p)\R = Pg{R) - (4r?sL + 3%g)^ - % ■ (78) 
r R R 



Here we made use of the velocity boundary condition, i.e. the first of Eqs. (|7^), and assumed the linear 
velocity profile inside the bubble which was introduced in Section 3, f = dr(j) = rR/R for r < R, 
which determines the gas viscosity term oc rjhG- As we pointed out already after Eq. (^), only the 
bulk viscosity of the gas contributes for a linear (homologous) velocity profile. 



Furthermore, as discussed in Section 4.1 , we employ the Navier-Stokes equation Eq. (53) as the 
underlying dynamical equation of motion of the gas and liquid fluids. Integrating over the radial 
coordinate r , we obtained the nonperturbative Eq. (|66|) before, which we presently rewrite as: 



dt^+-{dr^y] =I+h[PL]-^-d',{r(l))\R , (79) 
2 J B. pir 



where /i[-Pl] denotes the liquid enthalpy evaluated in Eq. (64) and we collected terms pertaining to 
the bubble interior into: 

I=-^RR + hG[PG] , (80) 



with the gas enthalpy ho from Eq. (|61|); note that ho = for a homogeneous density inside the 
bubble. We recall once again that the contribution of the bubble interior was evaluated assuming 
the linear velocity profile. The discussion following Eq. (^) furthermore implies that 1 = for the 
case of exactly homologous bubble interior dynamics (cf . Section 3 ) . Equations (^) , ( [78| ) , and ( |79|) 
completely determine the bubble motion, i.e. R{t) , as we shall demonstrate now. To begin with, we 
simplify the r.h.s. of Eq. (^) by using: 

fiLiPL] d^{r(j))\R = hL[PL-r]L-dr{r(j))\R] d^{r(j))\R „ , ^ +•.• , 81 

PLT r npLT Po + Pi 
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which follows from Eq. (|64|). Neglecting the small cross term (and higher order corrections) between 
viscosity and compressibility of the liquid as before, we eliminate Pl and replace the argument of 
hil- . .] with the help of the pressure boundary condition (^). The result is: 

hLlPciR) - i^VsL + 3vbG)iR/R) - i2a/R)]] = Hl , (82) 

which we henceforth abbreviate as indicated. 
In order to proceed, we calculate: 

dt4>\R = ifitret) + g'{tadv) " g' {tret)) r 

= ^{-cRR-ccP + 2g'{tadv))^ , (83) 

where we used the boundary condition ( |77[ ) in the second step. Inserting this result together with 
dr<P\R = R the above evaluation of back into Eq. (|79|), we obtain: 



cR + ^^\R-^g'{tadv)\R-^R^ =I + hL . (84) 

We have to eliminate cj), in order to arrive at the closed equation for R{t) which we are after. 

Taking the time derivative of Eq. (|8^), employing Eq. (^) once more, and solving the resulting 
equation for cp, one finds explicitly: 

c R ( i^, R. cR R^ 1 d , A g' g" , , 

where g' stands for g'(tadv)\R and g" similarly. Reinserting this result into Eq. (|83) yields: 



c 



■ R ■■ , R, cR, R, \ d , , \ 2 ,, 1 



from which equation we will derive our final results. 

First of all, we obtain a perturbative equation of motion for i?(t) by expanding the factor (l+i?/c)^^ 
in powers oi R/c and collecting terms. Taking up to second order corrections into account, the result 
is: 

RR[i-2- + 2-2 + 1 _ _ 2 _ ^ ^ ^/^ 

\ c c I I \ 6 c 6 c J c 

= (^l + |(l-|)|^(J + /.i) + 0((i?/c)3) . (87) 

It is straightforward to calculate higher order corrections in this case, if necessary. Equation (| 
when evaluated including only first order corrections in R/c, essentially agrees with our previous 
perturbative result, Eq. (|69|). To see this, note that dZ/dt generates a term ocR, using Eq. (|80|). The 



third order derivative can then be eliminated similarly as between Eqs. (p7D-(6£). 

However, a crucial difference remains in that presently we treat the external field exactly, which is 
represented here by the term oc g" , whereas in Section 4.1 the external field was treated perturbatively 
on the same footing as the sound field generated by the bubble wall motion. By comparison ofEqs. (|7D 
and (^5|) (the external field terms involving g ) we find here explicitly: 

--g"(t + R/c) = —sm(oj[t + R/c]) , 

C PL 
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for the example of an external standing plane wave pressure field discussed in Section 4.2 . We remark 
that R/c may vary systematically between less than about Ins and about 100ns during a nonlinear 
oscillation cycle with bubble parameters in the SBSL regime. Neglecting this retardation in the 
external field term, terms related to the bubble interior, and modifications due to the realistic liquid 
EoS and enthalpy, Eqs. ( p^ and (Bl), which enter our derivations, we again reproduce Eqs. (9) and 
(10) of Ref. i, if we truncate Eq. at 0{R/c) . 

Alternatively, instead of expanding Eq. (^6|), we simply multiply this equation by the denominator 
1 + R/c to obtain: 

RR{1 --) + U\l - - -(1 + -)g"{t + R/c) = {l + \l + hL) + -^il + hL) , (89) 
cz 6 c c c c c at 

which presents the nonperturbative bubble equation of motion; it agrees with the result of Ref. 
in the limit of the idealized adiabatic liquid EoS employed there and neglecting the bubble interior 
(except for the pressure boundary condition). 



5 Numerical Results 

We now turn to a numerical study of characteristic properties of a single air bubble in water at 
standard laboratory conditions (STP), i.e. at the ambient temperature T = 293 K and under the 
ambient pressure Pq = 1 atm = 1.01325 • 10^ kg/m • s^ . Solutions of the nonperturbative equation of 
motion (^) will serve as our standard, with which other results shall be compared. 

The following material 'constants' of water (STP) will be used H]: density pi = 10^ kg/m^ , surface 
tension a = 0.073 kg/s^ (water/air), sound velocity c = 1481 m/s . This latter value is about 2% higher 
than the one which follows from the adiabatic water EoS ( [6^ ) at STP with the parameters given in 
Ref. Q, n = 7 and Pi = 3 • 10^ kg/m • s'^ , which we employ in Eqs. (|6^ , (|8^). The shear and bulk 
viscosities of water (STP) are Q: r/^L = 1.002 • 10~^ kg/m • s , "q^L = 2.91r]sL . 

The 'standard' air bubble we study has an equilibrium radius Rq = 4- 10~^ m . The adiabatic index 
is 7 = 1.4 , characteristic of di-atomic molecules. The equilibrium radius and the particle number in 



the bubble, or equilibrium gas density po, are related through the vander Waals EoS (44) and pressure 



boundary condition ( [Tq ) with Pl = Pq (STP), which corresponds to the following set of parameters 
Q: A = 0.913 • 10^'^ , vander Waals excluded volume l^ra^ = 0.041/mole , ratio of equilibrium to hard 
core density po/pa = 0.0023 . For the air viscosities we take [||: rjhc = O.GtJsg = 1-08 • 10~^ kg/m • s , 
which indeed make a negligible contribution in the context of our study. 

The applied driving sound field, cf. Section 4.2, will be chosen to oscillate with the frequency 
u = Lo/I'K = 26 Khz in all examples. Its amplitude will be fixed at Pa = 1.35Po) except when stated 
otherwise. 

In the numerical calculations we consider the above specified parameters as prescribed constants, 
but keep in mind that e.g. the sound velocity and surface tension may change appreciably. The use 
of a constant adiabatic index has been questioned before since various transport effects leading 
to entropy production may become important. We follow Ref. [Q| in that we artificially increase all 
viscosities by a factor three, which has been found necessary to fit the standard Rayleigh-Plesset 
equation, i.e. R(t), to the experimental data obtained by lightscattering methods |jll|, p^ . 

In Fig. 1 we show the full cycle of the time dependent radius R{t) of the air bubble in water 
(full line), the bubble surface velocity R{t) (dashed line), and the sinusoidal external driving pressure 
Pa{^ = 0,t) according to Eq. ( [72[ ) (overlaid full line together with zero line, arbitrary units), which 
coincides with the relevant long- wavelength limit. We integrated numerically the nonperturbative 
Eq. (^) for a homologous bubble, in which case 2 = according to Eq. (|80|) . However, the gas 
pressure entering through Eq. (^) is evaluated here for a homogeneous bubble for the moment; we 
shall consider the effect of the inhomogeneity in Section 5.3 . We note that on the scale of Fig. 1 , 
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the R{t) curves calculated with the other equations derived in this work or the one used before in 
Refs. [Q, |l^, 15 1, for example, cannot be discriminated for our typical parameter set. Therefore, 
we will in the following pay more attention to R{t) , which is much more sensitive to the differences 
between the various equations of motion. 



5.1 Perturbative vs. Nonperturbative Equations of Motion 

The equations of motion which have been used to study driven gas bubbles in liquids in the SBSL 
parameter regime incorporated effects of the sound emission only in the lowest order in R/c so far, see 
e.g. Refs. lH, ^, |l^ and further references therein. In view of the high velocities of the bubble surface 
reached during the collapse phase, which may reach R/c ~ 1 (= Machl), as indicated in Fig. 1 , we 
study here how the various approximate treatments of sound damping compare to the nonperturbative 
Eq. (^). The bubble interior is treated here as explained in the context of Fig. 1 above. The individual 
sections in Fig. 2 show the surface velocity R{t) at the first bounce (left column) and the second bounce 
(second column) for different amplitudes of the driving pressure (top to bottom; Pa in units of Pq); 
the origin of the time-axis is arbitrarily chosen for the first bounce, the delay until the corresponding 
one for the second bounce is indicated in each case. The results are obtained from the nonperturbative 
equation (full line), and its lowest order in R/c perturbative expansion (dashed line), cf. Eq. (|87|), 
respectively. In the latter case also the retardation shifting the time argument of the driving pressure 
term, see Eq. (|88|), has been neglected, in conformity with earlier work. 

It is obvious from Fig. 2 that characteristic phase shifts result w.r.t. to the cycle timing set by the 
external driving pressure and between the first and subsequent bounce(s). Physically most important 
is the fact that within the perturbative treatment the maximal surface velocity during the collapse 
is overestimated by a considerable amount, e.g. ~ 250 ms~^ for the first bounce at Pa = 1.35Po • 
This implies that the collapse is significantly less violent than has been concluded from previous 
theoretical treatments of the bubble dynamics. In particular, quantitative SBSL estimates derived 
from models which couple the perturbatively corrected Rayleigh-Plesset equation, cf. Eq. (|6^), to a 



full hydrodynamic simulation of the bubble interior may need a revision |10| , 12, |I5|]. We are led to 
this conclusion also from the various other corrections to be discussed shortly. 

The origin of the retardation of the first bounce in the perturbative approach is shown in Fig. 3 
to be due to the neglect in the perturbative expansion of the time retardation effect. When the 



retardation is included as in Eq. (88), then the positions of the first perturbative and nonperturbative 
bounce coincide. However, the perturbative expansion including the retardation effect continues to 
introduce an overestimate of the collapse velocity. This is leading to a remaining phase shift of 15-20 
ns between first and subsequent bounce(s), as seen in the bottom portion of Fig. 3 . 



5.2 Homologous vs. Homogeneous Bubble Interior 

In order to demonstrate the influence of the behavior of the bubble interior on the overall dynamics 
of the coupled bubble-liquid system, beyond just providing the pressure resisting the collapse, we 
consider here how the homologous bubble description developed in Section 3 compares with the usual 
simpler homogeneous matter distribution model of the bubble interior. 

We exploit here the presence of a fixed point Eq. (^2|), where the density remains nearly constant. 
Presence of this fixed point cannot be expected in general, however, we presently make this assumption 
based on our numerical experience with the description of homologous and adiabatic changes of the 
interior distributions. We can evaluate Eq. ([45| ) conveniently integrating from ^^R to R. With P^: = 
P{p*), see Eq. (|^, which denotes the gas pressure of the corresponding homogeneous bubble, we 
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obtain the following approximate result for the homologous gas pressure at the bubble surface: 



Employing this expression for the gas pressure in the nonperturbative equation of motion (|89|), 
we find the results depicted in Fig. 4 (full line) , where also the behavior of a homogeneous bubble 
(dashed line) is shown for comparison. We observe only a rather small effect on R{t), except for a 
1.5 ns shift of the minimum. However, similarly to the effects illustrated in Fig. 2 , we find here another 
sizeable reduction of the maximal collapse speed (~ 200 ms~^). The maximum of the gas pressure at 
the homologous bubble surface is reduced by about 5% , as compared to the homogeneous one, which 
amounts to ~ 10^ atm . However, the maximal gas density at the surface is much less affected and 
reaches approximately one half liquid density in our calculational example. 

5.3 The Sound Field of a SBSL Bubble 

As a further application of our study of the sound field imposed on and rescattered by a typical 
SBSL bubble we evaluate it at a given distance away from the bubble, where it can be measured by 
a hydrophone, which preferably should be sensitive to the highest expected sound frequencies ||l], [Tl[ |. 

Perturbatively, it is most convenient to obtain the sound pressure using Eq. (|2^) or ( ^5|) in the 
form: 

P{r,t) = Pa{r,t) + - [p{R,i) - Pa{R,i)) , (91) 



where R = R{t) and i can be further evaluated with the help of Eq. (24) or Eq. (pq), whichever 
applies. The pressure difference on the r.h.s. here can then be computed easily from the solution of 
the perturbative equation of motion (|69|), cf. Section 5.1 , since it is implicitly part of it. 

The result at the lowest nontrivial order in R/c is shown in Fig. 5 , for r = 1 mm away from the 
center of the bubble. We clearly see the outgoing compression spikes riding on the sinusoidal driving 
pressure, which are caused by the bubble collapse and successive bounces. On the pressure scale of 
this figure (cutting off the first maximum) the nonperturbative result would be indistinguishable from 
the perturbative calculation. 

We calculate the pressure field nonperturbatively beginning with Eq. (|9l|), as before. However, 
we proceed in this case by using the pressure boundary condition, Eq. ([78|), in order to calculate 
P{R,t) = Pl{R). The term oc — r/L, in particular, can be evaluated in the long-wavelength limit 
(neglecting a cross term between viscosity and compressibility) to yield -t-2a;?7L-PA cos(a;t)/(3c^PL) , 
while all other terms are straightforward to obtain from the numerical solution of the nonperturbative 
Eq. (^); cf. the remarks explaining the calculation of Fig. 1 in Section 5.1 . 

In Fig. 6 we show the nonperturbative evaluation of the first sound pressure spike at r = 1 mm . 
For comparison, also the somewhat earlier and considerably stronger spike given by the perturbative 
calculation, cf. Fig. 5 , is shown. If we correct the calculated amplitude here for the geometric disper- 
sion, then the pressure at r = 0.6 /xm reaches 4.5-10^ atm at maximum. Applying the generic damping 
factor of 10~^ for the absorption of a (300 Mhz) pulse travelling 1 mm , we arrive approximately at the 
pressure amplitude measured experimentally for an SBSL bubble under similar conditions as assumed 
in our calculation [||, |ll|. However, another interesting aspect here is the rise time of the sound signal, 
i.e. from one-half to maximum amplitude, which takes only 40 ps (decay time 260 ps). This is about 
two orders of magnitude below what has been resolved in the above cited experiments. 

Finally, in view of the pressure spike results presented here, it seems worth while to check one of our 
basic assumptions, which is the linear acoustic approximation. Section 2.1. It is underlying the linear 
wave equation for the velocity potential, Eq. (|l3|), which forms the basis of our study of the emitted 
sound field. Whereas dispersion and absorption can be incorporated in a linearized approximation. 
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genuinely nonlinear effects are beyond our present scope The linear approximation requires that 
perturbations of the ambient state of the fluid are relatively small. Employing the adiabatic liquid 
EoS ( |63D for water and considering the maximum of the pressure amplitude obtained here, see Fig. 6 , 
we obtain a maximal compression and ratio of sound speeds: 

max-^ « 1.4 , max 4^ ~ 2.7 , (92) 
PL c[pl) 

respectively. Whereas the compression may appear to be quite tolerable, the temporary increase of 
the sound speed and its implications for the damping of the bubble motion, when the sound pulse is 
launched from the collapsing bubble surface, obviously deserve further study. 

Furthermore, one may wonder about the behavior of the water next to the bubble surface, when 
it is exposed to the "cold shock" indicated by our results, i.e. several 10^ atm pressure increase within 
less than 50 ps . 



6 Summary, Conclusions and Outlook 

The aim of our present work has been to reconsider the sound emission from the highly nonlinear, 
large amplitude motion of the interface between the gas inside and the liquid outside a cavitating 
bubble. Previous studies typically evaluated the radiated sound field in lowest order of a perturbation 
expansion in R/c, i.e. valid for slow bubble wall motion as compared to the sound velocity in the 
liquid. Generally, it is believed that R{t) reaches up to (or even exceeds) the sound velocity for 
externally driven bubbles in the parameter regime where sonoluminescence is observed experimentally 



This necessitates a nonperturbative treatment, such as ours, cf. Section 2.2 and see Eq. (28), 
in particular. We anticipate that the sound signal from a nonlinearly oscillating bubble wall may 
provide an important additional diagnostic tool. It should reflect the essential short-time scale(s) of 
this motion, which vary over several orders of magnitude 13], with current technology limiting the 
resolution in the (lOOMhz)"^ range, i.e. less than about 10 ns. This may be particularly valuable 
in cases where light scattering methods ^ do not truly reflect the motion of a sharply defined 
bubble wall or are not applicable at all, such as for liquid metals |p. 

We recall that based on Section 3 we describe the bubble interior using homologous profile functions 
for the density and velocity distributions. These allow for a more realistic study of the dynamics 
and particular properties of the high compression phase of strongly driven bubble oscillations, the 
importance of which has been shown earlier |lO| , |l2|| . Presently we employ only a simplified version 
of the semianalytic variational approach developed in Refs. |jl3|, In any case, our derivation of 
Eq. (1^) allows to incorporate easily any more precise model of the bubble interior; and it may be 
possible to extend this approach in order to account for the important effects of heat conductivity and 
mass diffusion. 

Our numerical examples obtained in Section 5 employ as another approximation the homologous 
solutions with the usual adiabatic vander Waals EoS. These EoS become questionable when the energy 
density reaches the ionization regime. The dynamic behavior of the gas mixtures within a rapidly 
oscillating bubble, here assumed to remain that of a di-atomic gas (air), deserves further study. 
In parallel to the cyclic in- and outgassing from the liquid into the bubble ("rectified diffusion"), 
it remains to be seen whether an essential amount of liquid vaporizes at the bubble surface and 
recondenses during an oscillation cycle. Related effects may help to explain the observed sensitivity of 
sonoluminescence to experimental parameters such as temperature in particular |l|, |l^]. In general, 
transport phenomena within the bubble and across the phase boundary have to be incorporated 
whenever the motion becomes fast compared to characteristic relaxation times Tj , ~ Aj/rj , where 
Aj denotes the scale of a corresponding gradient (e.g. mass density, partial pressure, temperature. 
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etc.). In particular, shock waves may be launched into the bubble interior 15, 16| or exterior |1C], 
the description of which is beyond the scope of this work. 

One may speculate whether or not the liquid (water) may be trapped in a metastable state w.r.t. 
solidification into a high-density phase (of ice) and in which form the corresponding binding energy 
would be released most efficiently. Presumably it stays at or close to the ambient temperature ||lO[| . 
Assuming an effective sound velocity of 2000 m/s, the 300 ps mean half width of the pressure pulse 
obtained in Section 5.3 corresponds to a spatial shell width of 0.6 /^m, i.e. about the vanderWaals 
hard core radius employed in our calculations. Near the collapsed state with i? ~ 0.6 /im such a 
shell contains about 10^^ water molecules, i.e. only about an order of magnitude more particles than 
assumed (in the gas) inside. Furthermore, the sound pulse amplitude decreases exponentially away 
from the bubble surface (besides geometric dispersion oc 1/r) due to absorption in the liquid Then, 
assuming an extremely conservative decay oc exp(— Ar [mm]/. 434) , one immediately estimates that 
the energy dissipated in this first shell (proportional to square of amplitude decrease) amounts to 
about 3/1000 of the pulse energy. This is more than an order of magnitude times the energy emitted 
in the form of visible light by a SBSL bubble Q . The small spatio-temporal extension of the region 
where this energy has to be dissipated, and how, definitely deserves further study. This should provide 
an improved starting point for the exploration of essential aspects of the violent bubble collapse which 
may help to elucidate the nature of sonoluminescence. We hope to come back to one or the other 
of these fascinating aspects of SBSL bubbles which are especially related to their sound field in our 
future work. 

We note that while in principle the effect of a dynamic treatment of the bubble interior on the 
behavior of R{t) , which we introduced in Section 5.2 and Fig. 4, should be part of a full hydrodynamic 



simulation |[T^, 12, |T^], the sound emission discussed in Section 5.1 had not been treated properly in any 
of the approaches prior to our work. In view of the numerical results presented in Section 5 it is obvious 
that the nonperturbative treatment of the emitted sound field as proposed here, which accounts for 
about 50% of the damping of the driven bubble oscillations, is mandatory. In particular, the standard 
perturbative calculations tend to overestimate the maximally reached surface velocity during the 
collapse phase considerably. Thus a crucial aspect of the "preparation phase" for the unknown light 
emission process can be described more accurately employing the nonperturbative equation of motion 
derived in Section 4.3 , Eq. (|89|) . This seems particularly relevant for detailed hydrodynamic studies 
of the bubble interior, in which the exterior has been described by the perturbative Rayleigh-Plesset 
equation before |l^, |12|, |15[. Our numerical results for the nonperturbative sound pulse emitted 
by an air bubble in water (Section 5.3) fit qualitatively the first experimental observations reported 
in Refs. |jl], 11|. On the other hand, as we pointed out at the end of Section 5.3, they also indicate 
the limitation of our present derivations. Whereas the bubble equation of motion, Eq. (|89|), is derived 



from the intrinsically nonlinear Navier-Stokes equation, in the form of Eq. (7£) of Section 4.3, our 
considerations of the emitted sound field are based on the acoustic approximation (see Section 2.2). 
Therefore, we still neglect important sound absorption and dispersion effects in the liquid. Using 
the realistic EoS for water, Eq. (^) in Section 4.1 , we indicated in Eqs. ( |92|) the compression of the 
liquid and the corresponding increase in the density dependent sound velocity which are induced by 
the outgoing pressure spike next to the bubble surface. It seems desirable to study in the future 
the importance of nonlinear effects on the propagation and fate of the extremely strong sound pulse 
emitted. 

There are published results on full hydrodynamic simulations of the coupled bubble-liquid system. 



which include especially the exterior fluid [16|. However, to the best of our knowledge, so far no 
attention has been paid to the unusual behavior of the liquid which may be caused by the pressure 
spike of 10^ . . . 10^ atm launched from the bubble surface, independently of whether there are shock 
waves generated in the bubble interior or not. We obtained such amplitudes in typical examples 
with generic SBSL parameters, where particularly their rise time of only about 40 ps seems quite 
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astounding, but also reminiscent of the shortness of the SBSL Ught pulses. 
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Figure Captions 

Fig. 1 : Radius R (full line) and surface velocity R (dashed line) as a function of time for an air bubble 
in water computed according to the nonperturbative Eq. (|89D. The time dependent driving pressure 
with Pa = 1.35Po is shown (overlaid together with zero line, arbitrary units). See main text for 
specification of other system parameters and further details. 

Fig. 2 : Surface velocity R{t) for the bubble of Fig. 1 at first bounce (left column) and second bounce 
(second column) for different amplitudes of driving pressure (top to bottom; in units of Pq); the 
delay between origins of the time-axes between the bounces is indicated. Results from the nonper- 
turbative Eq. (|89|) (full line) and the perturbative Eq. (|87|) (dashed line); see main text for further 
details. 

Fig. 3 : Same as Fig. 2 , however, employing the nonperturbative Eq. ( p9[ ) (full line) and perturbative 
Eq. (^), comprising time retardation (dashed line). 

Fig. 4 : Comparison of the radius R (top) and surface velocity R (bottom) for the bubble of Fig. 1 , 
using Eq. (|89|) , for a homologous (full line) and a homogeneous (dashed line) bubble. 

Fig. 5 : The pressure amplitude at r = 1mm from the center of the bubble of Fig. 1 , calculation 
including 0{R/c) corrections; see main text for further details. 

Fig. 6 : The first spike of the outgoing compression pulse launched by the collapsing bubble, see Fig. 5 . 
For comparison the result of the nonperturbative calculation (smaller, later pulse) is shown together 
with the perturbative one at 0{R/c) . 
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